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The dynamics of soliton pulses in the Nonlinear Schrodinger Equation (NLSE) driven by an 
external Traveling wave is studied analytically and numerically. The Hamiltonian structure of the 
system is used to show that, in the adiabatic approximation for a single soliton, the problem is 
integrable despite the large number of degrees of freedom. Fixed points of the system are found, 
and their linear stability is investigated. The fixed points correspond to a Doppler shifted resonance 
between the external wave and the soliton. The structure and topological changes of the phase 
space of the soliton parameters as functions of the strength of coupling are investigated. A physical 
derivation of the driven NLSE is given in the context of optical pulse propagation in asymmetric, 
twin-core optical fibers. The results can be applied to soliton stabilization and amplification. 
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I. INTRODUCTION 



Investigations of the externally driven Nonlinear Schrodinger (NLS) equation date back to the seminal work of 
Kaup and Newell Q| on the AC-driven damped NLS equation. That work was also one of the pioneering papers in 
which a perturbation method for solitons based on the Inverse Scattering Transform (1ST) technique was developed. 
Externally driven NLS equation arises in many applications, mainly in the context of solid state physics, such as long 
Josephson junctions M and charge density waves ||. The same equation describes plasmas driven by RF fields [Q. 
Much attention has been paid to the study of chaotic phenomena in the phase space of the soliton parameters of the 
driven system , and to the formation and stability of soliton states unique to the driven system Jr[j^] . Problems 
of generation of solitons via a coupling to external perturbations have been also investigated, for a homogeneous 
AC-drive ||] and traveling waves JhJ] . 

In the model system that we shall consider (see Eq. (|l|) below) the perturbation term does not introduce dissipation 
(the system remains conservative), and the uniform time-dependent driving is generalized to a, spatially dependent, 
driving field. The physical background of our model system is the nonlinear pulse propagation in optical fibers. In 
this case the driving force can assume, as will be shown in Sec. IV, the form of a Traveling wave coupled linearly 
to the NLS equation. The model system can be realized in the context of propagation of envelope optical pulses 
in twin-core optical fibers. The aim of this work is to study the effects of the driving field on the dynamics of the 
NLS-solitons. The emphasis is on weak couplings, and the main motivation is stabilizing and controlling solitons by 
an external traveling wave. 

We will limit ourselves to one-soliton pulses (that persist in an unperturbed NLS equation) and investigate the 
evolution of the parameters of the soliton when it is driven by a Traveling wave. This form of coupling preserves the 
Hamiltonian structure of the equation, thus making possible to employ Hamiltonian perturbation methods. We will 
show that, in the adiabatic approximation for a single soliton, the problem is integrable despite the large number 
of degrees of freedom. We will show that there exists a resonance, or phase locking, between the soliton and the 
driving field. We will see how this resonance is linked to the particle-like properties of solitons and how the resonance 
conditions are influenced by the strength of interaction. 

Transmission of solitons over long distances is essential for the use of solitons as digital bits in optical transmission 
lines . This has been achieved by use of Iridium doped fiber amplifiers . Here we will suggest an alternative 
scheme that employs a twin-core optical fiber. For identical fibers with pulses centered around the same central 
frequency, the twin-core system is also termed the Nonlinear Directional Coupler p3| , and this system has been 
extensively investigated (see |Q and references therein). Asymmetric twin-core optical fibers were also investigated. 
Numerical investigation of the asymmetric coupler as a means of performing logical gate operations with solitons was 
carried out in |jq| . In p6[ , "static" soliton states which can exist in these fibers with a limited asymmetry, were 
investigated in the model of two coupled NLS equations, where new types of solitons unique to the twin-core fiber 
were found. 

In our model, the two coupled NLS equations, which describe the evolution of optical pulses in asymmetric twin-core 
fibers, are reduced to a single NLS equation driven by an external traveling wave. We will see how the driving by 
an external wave can be used to stabilize soliton propagation in optical fibers, suggesting that this scheme can be 
used in transmission of solitons over long distances. We will also show how solitons can be amplified by the external 
Traveling wave with slowly varying parameters. 

The paper is organized as follows. In Sec. || we outline the Hamiltonian perturbation method for the NLS equation 
driven by an external traveling wave. In Sec. we show that the reduced problem is integrable and investigate a 
novel type of resonance (phase locking) between the soliton and the external wave. The linear stability analysis of this 
resonance is performed, and the structure and topological changes of the phase plane are investigated, analytically 



and numerically. A physical derivation of the driven NLS equation and possible applications are presented in Sec. IV. 
Sec. summarizes our results. 



II. HAMILTONIAN PERTURBATION METHOD FOR THE DRIVEN NLS EQUATION 

Consider the NLS equation coupled to an external traveling wave field: 

i^ + -^+2\ip\ 2 Tp = £eiq>[i(kx-ujt)} , (I) 

where we use dimensionless variables (see Sec. |Ey| for the corresponding physical units in the case of optical fibers). 
In Eq. (p]) e is the (normalized) strength of the coupling, and it is assumed to be small: e <C I. Also, ui and k are the 



2 



(normalized) frequency and wave number of the driving field, respectively. For specific physical models u> and k are 
related by an appropriate dispersion relation. 

The unperturbed (e = 0) version of Eq. ([l]) is completely integrable Jl7[ . and its most interesting solutions are 
solitons. The full one-soliton solution is given by: 
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where p, p, q and tp are free parameters. The quantity p defines the amplitude and width (1/p) of the soliton, p is 
the soliton's velocity, 2q/ p is the location of its "center of mass", and tp is its initial phase. 

Eq. (|l|) has Hamiltonian form, as it can be obtained from the variational derivative of a Hamiltonian: 
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where the Hamiltonian is given by 

H [ip] = 
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Here ip is the complex conjugate of ip, and and 9 are the real and imaginary parts of a complex number. 

For the unperturbed system, e = 0, the four parameters of the single-soliton solution (0) form a Hamiltonian 

dynamical system Jl8[ . For the variables q = q — ^t, and tp — tp — - — j — -t, the soliton amplitude p becomes the 
canonical momentum conjugate to the coordinate tp while the soliton velocity p is the canonical momentum conjugate 
to q. The (reduced) one-soliton Hamiltonian of the unperturbed system is given fl8| by: 
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Going back to the perturbed system, e ^ 0, we will employ the adiabatic approximation and investigate slow 
variations of the soliton parameters caused by the external field. In doing so we neglect radiation effects and possible 
formation of other solitons. This approach is physically motivated in the context of interaction of soliton pulses in 
optical fibers, where a typical initial condition is a single soliton solution of the unperturbed NLS equation. It is 
known that solitons are robust objects, especially in the case when the perturbations to the integrable system are 
Hamiltonian Therefore, for weak couplings, it is reasonable to expect that the main effect will be that of the 

solitons persisting, but slowly changing their parameters. 

It is advantageous to preserve the important Hamiltonian properties of the perturbed problem. Therefore, we will 
treat the perturbed system as a Hamiltonian system in the phase space of the one-soliton parameters driven by the 
external field. The driving term is obtained directly from the Hamiltonian (^J). Indeed, by inserting the unperturbed 
solution ip so \ with time- dependent parameters from Eq. (|^) into the Hamiltonian (Q), we obtain: 



H[ip sol \ = H (p,p, tp, q, t) = - | i>p- 
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The Hamiltonian (|6|) describes a dynamical system with two and a half degrees of freedom, as it is nontrivially 
coupled to an explicitly time-dependent driving force. As the dependence of the Hamiltonian upon tp and t enter 
only through a linear combination of these variables, a simple canonical transformation (see below) will eliminate 
the explicit time dependence, thus yielding an integral of motion: the new Hamiltonian. The new Hamiltonian with 
two degrees of freedom may still seem non-integrable, and one is tempted to look for chaos in this system. We will 
show, however, that because of the existence of an additional integral of motion, the one-soliton problem is actually 
integrable and can be fully investigated analytically. Formally, integrability occurs for any e. We expect, however, 
that the adiabatic approximation (that neglects radiation and possible creation of other solitons) will be valid only for 
small enough e. Therefore, wc will treat the driving term perturbatively and correspondingly require that e -C 1. On 
the other hand, one should assume that p > 1, to make the contribution of the nonlinear term in the NLS equation 
significant. This means that p ^> e which will be assumed to hold in the following. 
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III. PERTURBED ONE-SOLITON PHASE SPACE: INTEGRABILITY AND RESONANCE 



It can be easily checked that the perturbed one-soliton Hamiltonian (||) has an additional integral of motion, 

p(p - 2k) = const . (7) 
This integral follows from the complete, irreduced system, Eq. (pi), possessing an exact conservation law: 



ox 



const . 



(8) 



where = tp exp [—i [kx — ut)]. Indeed, inserting the one-soliton Ansatz (Q) for tfj in Eq. (^), one immediately obtains 
Eq -®- 

Note that the single-soliton solution form, Eq. (|7|) in our case, of exact integrals of the type (|8|), is only an 
approximation, valid as long as generation of radiation and formation of other solitons are ignored. Even so, it is a 
good approximation [ ^9|Jl4| when studying single soliton evolution under Hamiltonian perturbations. 

We now choose the integral of motion (Q) as a new momentum and make the corresponding canonical transformation 
in the one-soliton parameters' phase space. Simultaneously, we exploit the abovementioned symmetry in the time 
dependence of the one-soliton Hamiltonian and introduce the new phase $ (see below) . The generating function of 
this canonical transformation is: 



S{ Pl p,<$> 1 Q) = -{p-2k) P Q- p[<$> 

Therefore, the new canonical coordinates are given by 
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In the new coordinates the Hamiltonian becomes 



H(R,P,*,Q) = - 
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and the new momentum P is a constant of motion of the system. Therefore, the Hamiltonian ( |ll[ ) represents a system 
with one effective degree of freedom and is therefore integrable. This fact excludes any possibility of chaotic motion 
pot in the reduced, one-soliton system. 

Let us investigate possible resonances between the external wave and the soliton. In the language of the Hamiltonian 
(pi]), an exact resonance is related to a stable (elliptic) fixed point. Looking for fixed points, we should equate dH/dR 
and dH/d& both to zero. This yields two conditions. The first is 



$0 = 7T71 n = 0, 1 



(12) 



an exact condition, valid in all orders of e. Writing down the second condition, we will first limit ourselves to the zero 
order approximation with respect to s: 



R n = ± 



2(k 2 - u>) ± (4(fc 2 - u) 2 - P 2 ) ] 



1/2 



(13) 



The non-trivial, multi-valued resonance condition for the momentum R arises from the non-standard dependence of 
the Hamiltonian ( pi] ) on the momentum R. 

From Eq. © we obtain two conditions for the existence of a resonance between the driving field and soliton. First, 
k 2 > uj, which is a condition on the dispersion relation of the external wave alone. Second, 4 (fc 2 — un > P 2 , which 
is a condition on the parameters of the soliton and the external wave. 

The linear stability of the fixed points (|lj) and (13) is determined by the sign of the product FG, where G = 
d 2 H /dR 2 (R = Rq) (Ho is obtained from the Hamiltonian ( pl[ ) by setting e to 0), and F — e7rcosh _1 (itP/2Rq) 

The phase plane of the system ([II]) for a relatively weak coupling (e = 0.3) is shown in Fig. [|. The four-valued, 
in R, fixed points are clearly seen. The coordinates of the fixed points, found numerically, agree very well with the 
values given by Eqs. (12) and (PT 
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Fig. |2| shows the phase plane of the system ( |1 1\ ) for a stronger coupling, e = 0.6. One can see that a topology of the 
phase plane has changed. Four out of the eight fixed points have disappeared. In order to explain this bifurcation, 
we should modify the resonance condition, Eq. ( |l3| ) and take into account higher order corrections in s. 

In the first order in e the resonance condition for R becomes: 

47r 2 Psinh(^ N ) 

R i -A(k 2 -u J )R 2 +P 2 T e ¥*J-=Q, (14) 

Rn cosh 2 f ^ 



where the values for Rq should be taken form the zeroth approximation (O) . From Eq. ( |l4| ) we can see that indeed 
not all the fixed points which were present when e was small (see Eqs. (|12|) and (|l3|)) still persist when the coupling 
is increased. For a given P, when 



PR cosh 2 (§^) 
e > 



4-7T 2 sinh 



two of the resonant values of R do not exist anymore, thus leading to a topological change in the phase plane, as 
observed in Fig. 0. 

Now let us discuss the resonance condition in physical terms. Transforming back to the "lab" coordinates, using 
Eq. (|Io|), we can write the resonance condition $ = as 

. „_„. (15) 



dt V R 

Using the definition of P in Eq. (fLOj), Eq. (|15|) can be rewritten as 



dt \R 2 

But P is a constant of motion of our reduced system, so P = 0. Also, R — at resonance. Therefore, Eq. © 
reduces to 

¥>- (Sr^)^ = - (17) 

Now, g = Rv s /2, where v s in the center-of-mass velocity of the envelope of the soliton. Therefore, defining k = p/2 as 
the soliton internal wave wave number (see Eq. (0)), and uj c = ip (in the leading order in e), the resonance condition 
becomes 



U) c — KV S = LU — kv s . (18) 

From Eq. (|l8|) we can see that the resonance described by the reduced Hamiltonian (|ll]) is a Doppler shifted 
resonance between two waves: the external (pumping) wave with the wave number k and frequency to, and the carrier 
wave of the soliton with the wave number k and frequency lo c . It is interesting that there are two Doppler shifts in 
the resonance condition. The first Doppler shift enters the right hand side of Eq. (|l8|), and it is by the center-of- 
mass velocity of the soliton. The second Doppler shift, entering the left hand side, is less intuitive, and is also by 
the center-of-mass velocity of the soliton, but this time with the soliton carrier wave's wave number. This problem 
provides us with another non-trivial example of particle-like properties of solitons. 

Let us look at the limiting case of fc — > 0, in which the system should reduce to well known resonance with a 
homogeneous AC-drive jij. For oj < and p — (which are the parameter values studied in |ij), the constant of 
motion P = and the resonant values of the soliton's amplitude are reduced, see Eq. (|l3|), to p = R = ±\/2oj, as 
obtained in jij. For uj > 0, we find from Eq. ( p"3| ) that since k 2 > u> must hold there will be no resonance for k — 0. 
So that in this case the resonance is unique to a coupling to an external travelling wave. 

Another interesting limiting case is when the values of the wave vectors of the driving field and the internal soliton 
wave are close to each other. That is k ~ n = p/2. This translates (see the transformation (fio|)) to P — > 0. From Eq. 
( |l4| ) we can see that in this limiting case the reduced system will be exactly at the bifurcation point, for all values of 
e. Note that for k = 0, also v s — so that the resonance (ttq) is reduced to a simple resonance. 
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IV. PHYSICAL MODEL AND APPLICATIONS 



In this Section we will present a physical derivation of Eq. (|lj). We will start with the equations for the envel opes o f 
pulses in the cores of two adjoining, closely spaced, non-identical, single-mode fibers (twin-core optical fibers) 



+ "^Jr + 2 l - ^ 1 ! 2 ^1 + ^"12^2 exp —i (kx - LOtj 



0, 



(19a) 



dj>2 _ a 

dx dt 



fa^jr +21^1^2 + ^1 exp 



i [kx — u)t\ 



= 0. 



(19b) 



The coordinates x and t in Eqs. ( jig) ) are written in the "soliton units" p3] corresponding to Eq. ( |l9a| ). The couplings 
ai2 and a2i result from the overlapping of the evanescent fields of the transverse fiber modes with the fields in the 
adjoining fiber cores. Since the fibers are not identical, the coupling is not symmetric, i.e. an ^ a-ix- It is assumed 
that the transverse fiber mode is not affected by the proximity of the adjoining fiber, and by the identical transverse 
mode in it. We also assume that the interaction term arising from the cross phase modulation (term proportional to 

|V>i| 2 ">p3-i, i — 1,2) can be neglected. 5 — ("fi/^) 1 ^ 2 is the ratio of the nonlinearity strengths in the two fibers, where 



It = 



cA 



off 



i = 1,2. 



(20) 



Af s is the effective core area (which scales like pf, pi being the fiber core radius), n 2 is the Kerr coefficient, c is the 
speed of light, and uii is the carrier frequency in each fiber. The amplitudes ipi of the pulses are scaled, following [p3| , 
to (7i//3) 1 ^ 2 To, where [3 is the dispersion coefficient of the pulses in Eq. (19a), and Tq is the pulse width. 

It follows from Eq. (|2^) that the inequality 8^1 may result from the fibers having different radii, in which case S is 
the ratio between the radii of the fiber cores of the two fibers. Also, if the fibers are centered around different central 
frequencies u>i, then 5 = (lui/^) 1 ^ 2 ■ Let us continue our discussion of the different coefficien ts in Eq. jis|). The 



coefficient /3i is a measure of the difference in the group velocity in Eq. (19t) from that in Eq. (19a). The coefficient 
02 is the ratio of the dispersion coefficients of the two fibers, ft ^ 1 may result from the fibers not having the same 
transverse wave numbers, or from the pulses in the fibers being centered around different central frequencies. The 
different carrier frequencies and/or transverse wave numbers also lead to the fibers having different phase velocities. 
This fact results in the oscillatory term in the interaction, with k and u) being the mismatches in the wave-number 
and frequency, respectively. This oscillatory term arises in the evaluation of the overlap integral of the transverse 
modes of the two fibers. Notice that, as the coupling in Eqs. (19) is asymmetric, it is not possible in general to cast 
the system in a Hamiltonian form as it was done in p^] . We will show below, however, that in a certain limit the 
system (|l9|) can be reduced (see Eq. (p2])) to an equation possessing Hamiltonian structure. 

W e assume that the interaction term (the r.h.s.) in Eq. ( 19a ) is much larger than the interaction term in Eq. 
( |19b| ): 8a\2 S> a2i/8. This condition implies that 6 S> 1 which occurs when the ratio between the radii of the two 
fibers is large, and when the carrier frequencies are not the same. Under these conditions the interaction term in Eq. 
( |19b| ) c an b e neglected, and Eq. (19b) is decoupled from Eq. (19a) in the sense that it only enters as a d rivin g term 
in Eq. (19a), while there is no back action. Now, if we further assume that the pulses descri bed b y Eq. ( |l9b| ) are in 
the positive (normal) dispersion regime then there is no modulational instability [ p4| in Eq. (19b), a nd s table linear 
dispersive waves can propagate in the fiber. We are interested in the small amplitude limit of Eq. ( |l9b[ ), when the 
pulses are jus t linear waves. In this case we can drop the term arising from the Kerr nonlinearity. Then the set of 
equations (19) can be written as: 



l ~Q^r + -QtfT + 2 IW "01 +£"02 exp 



—i [kx — cut 
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(21a) 



dip2 
dx 



chp2 
dt 



■A 



dt 2 



(21b) 



where e = 5a\2. Therefore t he sy stem ( |l9[ ) reduces to a single NLS equation driven by an external Traveling wave. 
The equation for Eq. (21a), can be written (omitting the indices) as: 



G 



i^ + ^- + 2\^\ 2 iP + £exp[i(Kx-{lt)} = 0, (22) 



where the dispersion relation f2 (K) is given by the equation for ip2 hi Eq. ( piq ) together with the phase mismatch k 
and cD from Eq. (|l9|). 

Eq. ( p2] ) is equivalent to Eq. ([!]) with a; and t interchanged, K = —lu, and fl = —k. Therefore we can apply the 



results obtained in Sec. Ill to the analysis of the dynamics of soliton pulses in twin-core fibers, under the conditions 
that lead to Eq. (||). 

One application of our results is soliton phase locking. The resonance conditions [see Eqs. ( |l2| ) and fllq)] correspond 
to the soliton's parameters being phased locked to the driving field. For the stable fixed points this allows fixing of 
the soliton parameters (see also Fig. |l|). Specifically, the soliton amplitude p, see Eq. (||),_is constant. The values of 



the soliton parameters for which the resonance condition is satisfied are given, using Eq. (JlOj) , by Eqs. (f2) and (|13|), 
and by the value of the constant of motion, Eq. (Q), which is determined by the initial conditions of the soliton pulse 
and by the parameters of the driving field. Furthermore, the resonance is sustained, as can be seen from Fig. |^, also 
for large values of coupling (e in Eq. (|l|)), only that the number of resonances is decreased (see Eq. (14)). 



The phase locking can be used in order to stabilize solitons. In any real soliton transmission system there exists 
dissipation due to fiber losses. The dissipation will result in a decrease of the solitons' amplitude. The dissipation can 
be incorporated into the NLS equation by the addition p3 of a perturbation term in the form —iTip. Although the 
dissipative term cannot be directly included in our Hamiltonian perturbation approach, we can anticipate, in analogy 
to the problem of a driven damped oscillator, that by keeping the soliton's parameters in resonance with the driving 
field we will be able to overcome the effect of the losses in the system. In this case, the elliptic fixed points of Figs, [j] 
and H will become attracting points. 

Another application for which the coupled system can be used is soliton amplification. The amplitude (/? in Eq. 
(^J)) of the soliton can be increased for pulses which parameters correspond to periodic orbits surrounding stable fixed 
points (see Figs, [j] and |J). The period of oscillations for the stable orbits is long and, for small oscillations, is given by 
e l / 2 FG, where F and G were defined following Eq. (|l3|). By a proper choice of the initial conditions and interaction 
length the soliton amplitude can be increased by performing a half period of nonlinear oscillation around a stable 
fixed point. 

Another mechanism, by which the amplitude of solitons can be increased more significantly, is the "Dynamic 
Autoresonance" ||[L(^|5). The frequency of nonlinear oscillations depends on their amplitude. Therefore, the infinite 
growth of the amplitude, obtained for a linear, dissipation free oscillator driven by a resonant external force, is not 
possible for a nonlinear oscillator with constant parameters. By "chirping" adiabatically the driving field's frequency, 
one can preserve the phase plane structure so that the oscillator will continue to perform nonlinear oscillations around 
the (time dependent) resonant value of the action variable, R in Eq. (^) in our case. The soliton parameters will 
remain phase locked in resonance with the (slowly varying) driving field. Since R is the soliton amplitude, see Eq. 
(O), this mechanism provides a means for substantially increasing the soliton amplitude in a resonant manner. 



V. SUMMARY 



We have investigated the evolution of single-soliton pulses of the Nonlinear Schrodinger (NLS) equation driven by 
an external traveling wave field. This system, even though not integrable, is still Hamiltonian. Using the Hamiltonian 
structure and adiabatic approximation for a single soliton, we reduced the perturbed NLS equation to a two-and-a- 
half dimensional Hamiltonian system in the phase space of the parameters of the single soliton solution. One integral 
of motion in this system results from the fact that the time-dependence drops out in a rotating reference frame. 
An additional integral of motion is a consequence of an exact integral of motion in the complete, unreduced partial 
differential equation. Therefore, the reduced system becomes effectively one-dimensional and therefore integrable. 
Physically, the reduced system represents a nonlinear oscillator (with an unusual form of the Hamiltonian) coupled 
to an external harmonic force. The phase plane of this system was investigated analytically and numerically, and a 
good agreement between these two was found. As the coupling strength increases, there occurs a bifurcation in the 
phase plane of the reduced system. This bifurcation has been explained analytically. 

We gave a physical motivation to this model by showing that dynamics of soliton pulses in twin-core, non-identical, 
single-mode optical fibers can be reduced to a system in which the dynamics in one core are governed by a NLS 
equation driven by a linear Traveling wave propagating in the adjoining core. In this regime, one can neglect the back 
action of the nonlinear wave on the linear wave. 

Finally, we discussed possible applications of our results to stabilization and amplification of soliton pulses in the 
asymmetric, twin-core optical fiber. 
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FIG. 1. The phase plane of the Hamiltonian system ([Tl]) for a weak coupling. The parameters are e — 0.3, P = 3.0, k = 2.0 
and ui = 2.0 . The resonances corresponding to the four-valued solution of Eq. ( |l3| ) are clearly seen. Note that the stability of 
the resonant points is interchanged for positive and negative values of R. 



FIG. 2. The phase plane of the Hamiltonian system, Eq. (|ll|), f° r a stronger coupling, e = 0.6, The parameters are the same 
as in Fig. [I], except that e = 0.6. A bifurcation has occurred. In accordance with Eq. (|l4|), for each of the resonant values of 
R there are no longer two resonant phases, one for a stable (elliptic) point and one for an unstable (hyperbolic) point. Only 
one fixed point is left for each value of R. The phase plane topology has changed accordingly. 
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